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We  report  the  results  of  a  los^energy  electroij&iiffraction  analysis  of  th^atomic  structure  of  the 
{ 3 10 1  surface  of  body-centered -cubic  iron.  Damped  relaxations  from  bulklike  structure  are  found  in 
the  directions  both  parallel  and  perpendicular  to  the  surface  and  extending  fouf layers  deep  into  the 
lattice.  The  surface  structure  is  as  follows:  firstiptertoyer  spacing  rf  ”= 0.76ffD.0^A  (bulk  value jp 
0.906  A),  second  interlayer  spacing  d^  =  1.02&^.03  A,  third  intedayer  spacing  dU= 0.87,^0.04  A, 
change  of  first-to-sec^id  layer  registry  shiftjSo*  =0. 13^0.05  A,  and  change  of  second-to-third 
layer  registry  shift  (4^23  =0.03(^6.05  A.  Both  changes  of  registry  shifts  are  in  the  same  direction, 
along  a  (310)  axis,  and  cause  a  decrease  in  the  difference  between  nearest-  and  next-nearest- 
neighbor  distances  between  atoms  in  adjacent  layers.  £,  1 
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I.  INTRODUCTION  (It  ' 

Renewed  interest  in  surface  structure  of  metals  has 
developed  in  the  last  few  years  as  advances  in  both  surface 
theory  and  experimental  techniques  have  revealed  consid¬ 
erable  complexity  in  the  relaxation  of  the  outermost  atom¬ 
ic  layers.  Calculations  have  predicted,1'2  and  a  number  of 
experiments  have  shown,3-10  deviations  from  bulklike 
structure  extending  several  layers  deep.  Structures  for 
low-index  faces  exhibit  a  characteristic  damped  oscillatory 
relaxation  of  interlayer  spacings  (i.e.,  changing  sign  from 
layer  to  layer).  Recently,  with  the  case  of  Fe{21 1  ),10  we 
have  demonstrated  that  on  surfaces  less  symmetrical  than 
the  commonly  studied  (001  j,  { 1 10} ,  and  (111)  surfaces, 
multilayer  relaxations  involving  motion  of  atoms  parallel 
to  the  surface  occur  and,  for  the  (211)  surface  at  least, 
these  relaxations  are  also  oscillatory.  However,  the  pat¬ 
tern  of  relaxations  on  surfaces  of  lower  symmetry,  i.e.,  the 
signs  of  the  interlayer  relaxations  and  the  registry  shifts, 
appear  to  vary  from  surface  to  surface  for  the  three  cases 
which  have  now  been  studied:  Fe(211j,  Fe(210j,11  and 
Fe(310).  It  is  this  complexity,  which  is  a  function  of 
both  electronic  and  atomic  structure  near  the  surface  and 
which  provides  a  test  for  theories  of  metal  surfaces,  that 
makes  these  cases  very  interesting.  By  essentially  dou¬ 
bling  the  amount  of  experimental  information — the  relax¬ 
ation  of  each  atomic  layer  now  has  two  components,  one 
parallel  and  one  perpendicular  to  the  surface — the  test  of 
theories  of  surface  structure  by  comparison  of  observed 
and  predicted  structures  is  made  more  stringent.  Fe{  310} 
was  chosen  for  study  because  it  is  a  quite  open,  or  rough, 
surface  (packing  fraction12  0.3725)  with  low  symmetry 
(only  one  mirror  plane);  these  properties  make  it  a  favor¬ 
able  candidate  for  multilayer  relaxation.  Section  II  de¬ 
scribes  the  experimental  methods,  Sec.  Ill  briefly  outlines 
the  calculations,  and  Sec.  IV  presents  the  structure 
analysis  of  the  low-energy  electron-diffraction  (LEED)  in¬ 
tensity  data  by  means  of  experiment-theory  comparison, 
using  the  reliability  factor  (r  factor)  of  Zanazzi  and 
Jona.13 


II.  EXPERIMENTAL  SETUP 

The  Fe{  3 10}  sample  was  spark-cut  from  a  single¬ 
crystal  specimen  grown  by  strain-annealing  of  ultrapure 
source  material  provided  by  the  American  Iron  and  Steel 
Institute  (for  details  see  Ref.  14).  Laue  photographs  were 
used  to  orient  the  sample  to  within  ±0.5“  of  the  (310J 
plane.  The  crystal  was  then  ground  and  mechanically  pol¬ 
ished  with  water-based  alumina  slurries  (smallest  bead 
size  0.3  fxm). 

The  following  steps  were  required  to  obtain  a  clean  sur¬ 
face  in  ultrahigh  vacuum  (base  pressure  about  7X10“ 11 
Torr):  (i)  room -temperature  Ar+  bombardment  (2—4  ft  A 
at  400  eV)  for  2  h,  (ii)  Ar+  bombardment  .vith  sample 
temperature  raised  to  850  “C  for  2  h,  followed  by  a  1-h  an¬ 
neal  at  the  same  temperature,  and  (iii)  a  series  of  Ar+ 
bombardments  for  2—3  h  each  with  sample  temperature 
at  400— 500  “C  followed  by  1-h  anneal  at  600— 650  “C  (to¬ 
tal  time  of  argon  treatments  approximately  25  h).  Auger 
spectroscopy  was  used  to  monitor  surface  cleanness; 
sulfur,  oxygen,  and  carbon  were  the  main  impurities  be¬ 
fore  cleaning.  After  cleaning,  surface  cleanness  was  com¬ 
parable  to  our  previous  iron  work.15 

Thirty-three  low-energy  electron-diffraction  (LEED) 
spectra  (intensity-versus-energy  curves,  21  nondegenerate) 
were  collected  by  means  of  a  spot-photometer,  the  sample 
being  held  at  room  temperature  for  all  measurements: 
fourteen  spectra  at  normal  incidence  (10,  01,  10,  01,  11, 
1 1,  21,  12,  22,  22,  32,  23,  32,  and  23),  ten  spectra  at  0=7°, 
4>=-  107.5“  (00,  10,  01,  10,  01,  Tl,  12,  21,  23,  and  32), 
and  nine  spectra  at  0=14“,  <f>=  — 107.5“  (00,  10,  01,  10, 
01,  11,  21,  12,  and  23).  Angles  and  beam  indices  follow 
the  convention  of  Zanazzi  et  a/.1 6  Figure  1  shows  a 
schematic  representation  of  the  LEED  pattern. 

The  method  used  to  determine  the  condition  of  normal 
incidence  of  the  incoming  electron  beam  on  the  sample 
was  the  following:  (i)  orient  the  crystal  so  that  the  incom¬ 
ing  beam  is  in  the  ( 100)  mirror  plane  normal  to  the  sur¬ 
face  by  adjusting  the  sample  position  until  the  beams  with 
indives  xy  and  yx  in  Fig.  1  (with  x,y  equal  to  0, 1,2,3,. . .) 
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FIG.  1.  Schematic  LEED  pattern  from  Fe(310)  for  0^tO* 
and  —  107.5*. 


give,  as  nearly  as  possible,  identical  /-  V  spectra,  and  (ii) 
keeping  the  incoming  beam  in  the  { 100)  mirror  plane,  ro¬ 
tate  the  sample  by  various  angles  6  such  that  the  00  scat¬ 
tered  beam  appears  on  the  fluorescent  screen  on  both  sides 
of  the  electron  gun  (located  at  the  center  of  the  screen) 
and  record  the  corresponding  I-V  spectra  of  the  00  beam 
at  each  angle;  by  symmetry,  for  equal  6  on  opposite  sides 
of  the  LEED  gun  the  00  spectra  will  be  identical  and  can 
thus  be  used  to  determine  the  condition  0=0*. 


HI.  CALCULATIONS 

Dynamical  LEED  intensity  calculations  were  done  with 
the  computer  program  change  (described  in  Refs.  17  and 
18).  The  Fe  potential  was  the  same  as  used  successfully 
for  studies  of  other  iron  surfaces  (for  example,  Legg 
et  a/.19).  The  inner  potential  was  taken  as  energy  in¬ 
dependent;  the  imaginary  part  P  of  the  potential  was  set 
equal  to  4  eV  and  the  real  part  V0  was  to  be  determined 
from  the  intensity  analysis;  the  initial  value  was  fixed  at 
—  11.5  eV.  The  mean  atomic  vibrational  amplitude  was 
taken  as  ((u2))1/2=0.115  A  and  eight  phase  shifts  were 


used  to  represent  the  ion  core  scattering.  Occasionally, 
some  troublesome  numerical  instabilities  were  encoun¬ 
tered  (blowups,  i.e.,  very  large  intensity  values,  and  un¬ 
physical  kinks  in  the  calculated  curves),  presumably  due 
to  the  small  interlayer  spacings  of  Fe{  3 10)  (dbuik  =0.906 
A).  The  main  remedy  was  to  bunch  together  layers  into 
slabs  three  or  four  layers  thick,  with  scattering  within  the 
slabs  done  in  the  angular  momentum  basis  (which  has  no 
limitations  in  handling  closely  spaced  layers)  and  scatter¬ 
ing  between  slabs  done  in  the  plane-wave  basis.  Up  to  65 
plane  waves,  or  beams,  were  used  in  the  calculations.  At 
245  eV,  65  beams  are  propagating. 

IV.  STRUCTURE  ANALYSIS 

Top  and  side  views  of  the  outer  atomic  layers  of  a  bulk¬ 
like  Fej  310)  surface  are  shown  in  Fig.  2,  and  various  sur¬ 
face  parameters  are  defined  in  Table  I.  The  surface  is 
quite  open,  with  a  packing  fraction12  of  0.3725  (the  close- 
packed  bcc  1 110)  face  has  a  packing  fraction  of  0.8330). 
The  surface  unit  mesh  is  a  parallelogram  formed  by  mesh 
vectors  a2,  directed  along  a  <  100)  axis  and  of  length  2.87 
A,  and  aj,  directed  along  a  (311)  axis  and  of  length  4.75 
A.  The  interlayer  translation  vector  c  is  directed  along  a 
close-packed  (111)  axis  and  the  bulk  interlayer  spacing  is 
0.906  A.  The  bulklike  registry  of  one  layer  with  respect 
to  the  next-deepest  layer  is  given  by  the  relation 
cp«=T®i  +  ^a2=TaV^0^+T‘1T.  where  c^,  is  the 
projection  of  the  interlayer  translation  vector  c  onto  the 
(310)  plane,  and  x  and  y  are  defined  in  Table  I.  The  sur¬ 
face  has  one  mirror  plane,  the  j  100)  plane,  but  there  are 
no  rotational  symmetries. 

The  low  symmetry  of  the  (310)  surface  (see  Fig.  3) 
leads  to  the  expectation  that  the  surface  relaxation  will  ex¬ 
hibit  registry  shifts,  i.e.,  parallel  relaxation,  in  addition  to 
the  usual  changes  in  interlayer  spacings,  as  has  been  found 
for  Fe(211).10  For  example,  if  atoms  of  the  first  atomic 
layer  were  to  move  in  the  direction  opposite  to  the  axis  la¬ 
beled  (310)  in  the  top  view  of  Fig.  2  by  a  distance  of  0.23 


SIDE  VIEW  {310} 


FIG.  2.  Top  and  side  views  of  the  undistorted  Fef  310)  surface.  Circles  of  equal  thickness  represent  coplanar  atoms;  for  the  top 
view  circles  of  decreasing  thickness  indicate  progressively  deeper  layers.  All  distances  are  in  A. 
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TABLE  I.  Surface  lattice  vectors  at  and  a2,  parallel  com¬ 
ponent  of  the  interlayer  translation  vector,  interlayer  dis¬ 
tance  d,  and  mirror  plane  m  for  the  bulklike  Fe|310)  surface. 
$  is  along  (310)  in  Fig.  2  and  perpendicular  to  p;  $  is  along  a2 
and  parallel  to  ( 100)  in  Fig.  2;  z  is  directed  along  inward  sur¬ 
face  normal,  a  =2.  S66  A. 
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A,  then  each  top-layer  atom  would  have  three  nearest 
neighbors  in  the  second  layer  as  compared  to  two  for  the 
bulklike  structure.  As  a  result  of  the  shift,  however,  the 
nearest-neighbor  distance  between  atoms  in  the  first  and 
third  layers  would  then  be  decreased  and  the  correspond- 


FIG.  3.  Schematic  representation  of  the  first  four  layers  of 
the  Fe(310)  surface.  Solid  circles:  first  layer.  Triangles: 
second  layer.  Squares:  third  layer.  Hexagons:  fourth  layer. 


TABLE  II.  Results  of  calculations  for  6=0*  of  various  structural  models:  di2,  d23,  and  are  the 
first  three  interlayer  spacings,  a  i2  and  a23  are  the  first  two  interlayer  registry  parameters  (see  text).  The 
agreement  between  theory  and  experiment  is  given  by  the  reliability  factor  r.  All  distances  are  in  A. 


Run 

dn 

d2j 

d» 

°I2 

«2J 

r 

hbulk) 

0.906 

0.906 

0.906 

1.813 

1.813 

0.292 

2 

0.880 

0.906 

0.906 

1.813 

1.813 

0.280 

3 

0.854 

0.906 

0.906 

1.813 

1.813 

0.278 

4 

0.854 

0.933 

0.906 

1.813 

1.813 

0.219 

5 

0.854 

0.933 

0.906 

2.015 

1.813 

0.199 

6 

0.854 

0.933 

0.906 

1.611 

1.813 

0.250 

7 

0.854 

0.933 

0.906 

2.015 

1.914 

0.191 

8 

0.854 

0.933 

0.880 

2.015 

1.914 

0.171 

9 

0.822 

0.933 

0.880 

2.015 

1.914 

0.152 

10 

0.822 

0.933 

0.922 

2.015 

1.914 

0.162 

11 

0.822 

0.933 

0.943 

2.015 

1.914 

0.151 

12 

0.801 

0.933 

0.943 

2.015 

1.914 

0.137 

13 

0.801 

0.933 

0.922 

2.015 

1.914 

0.152 

14 

0.801 

0.933 

0.880 

2.015 

1.914 

0.144 

15 

0.779 

0.933 

0.943 

2.015 

1.914 

0.130 

16 

0.811 

0.943 

0.943 

2.015 

1.914 

0.144 

17 

0.779 

0.949 

0.938 

2.015 

1.914 

0.127 

18 

0.779 

0.965 

0.938 

2.015 

1.914 

0.129 

19 

0.779 

0.949 

0.880 

2.015 

1.914 

0.136 

20 

0.779 

0.949 

0.880 

1.964 

1.914 

0.125 

21 

0.779 

0.949 

0.943 

1.914 

1.813 

0.136 

22 

0.779 

0.949 

0.901 

1.914 

1.813 

0.124 

23 

0.779 

0.949 

0.880 

1.914 

1.813 

0.122 

24 

0.779 

0.949 

0.880 

1.939 

1.813 

0.120 

25 

0.779 

0.949 

0.880 

1.889 

1.813 

0.125 

26 

0.779 

0.949 

0.880 

1.864 

1.813 

0.130 

27 

0.759 

1.U14 

0.880 

1.994 

1.813 

0.102 

28 

0.759 

1.014 

0.880 

1.964 

1.813 

0.100 

29 

0.759 

1.014 

0.880 

1.934 

1.813 

0.100 

30 

0.759 

1.014 

0.880 

1.903 

1.813 

0.102 

31 

0.759 

1.014 

0.880 

1.949 

1.934 

0.113 

32 

0.759 

1.014 

0.880 

1.949 

1.874 

0.100 

33 

0.759 

1.014 

0.880 

1.949 

1.813 

0.100 

34 

0.759 

1.014 

0.880 

1.949 

1.752 

0.113 
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FIG.  4.  Variation  of  the  r  factor  for  thedata  set  at  6=0“  as  a 
function  of  dn  and  d23  with  =0.880  A,  0|2  =  1.934  A,  and 
V0=  — 10.5  eV. 


ing  next-nearest-neighbor  distances  increased,  which 
would  be  expected  to  limit  the  top-layer  registry  shift. 
For  parallel  shifts  of  deeper  layers  it  is  unclear  even  what 
the  signs  of  the  shifts  might  be.  The  purpose  of  the 
LEED  intensity  analysis  carried  out  in  this  work  was  to 
determine  the  signs  and  magnitudes  of  the  registry  shifts 
and  the  relaxations  of  the  interlayer  spacings. 

The  structural  parameters  that  were  varied  in  the  calcu¬ 
lations  were  the  first  three  interlayer  spacings,  di2,  d23, 
and  d M,  and  the  first  two  registry  shifts  au  and  a2j. 
(The  registry  shift  between  successive  layers  i  and  j  is  de¬ 
fined  by  ~cfmTij  =  \l2+aij£,  where  is  the  projection 
of  the  interlayer  translation  vector  onto  the  {310}  plane. 
See  Table  I  and  Fig.  2  for  definitions.  The  bulk  value  of 
atj  is  1.813  A.)  Beginning  with  the  bulklike  structure,  all 
di  equal  to  0.906  A  and  all  ai}  equal  to  1.813  A,  a  series 
of  calculations  for  normal  incidence  were  performed 


U  Bl!  (I  60  (BULK) 
184  1.94  EXPT 

aO  ®S 


(BULK) 

(0.906) 

(0.906) 
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FIG.  5.  Schematic  representation  with  optimized  surface 
parameters  for  Fe| 310 j .  All  distances  are  in  A. 


which  involved  independent  and  relatively  large  variations 
of  the  five  structural  parameters.  A  sampling  of  the  cal¬ 
culations  is  shown  in  Table  II,  labeled  with  run  numbers 
1—23.  The  models  that  gave  the  best  agreement  with  ex¬ 
periment  indicated  that  d12,  d2 3,  and  di4  were  contracted 
(  —  15%),  expanded  (—5%),  and  contracted  (—3%), 
respectively,  compared  to  bulklike  values,  that  a12  in¬ 
creased  by  about  5%  over  the  bulk  value  and  a2 3  was 
nearly  unchanged. 

Next,  with  di2,  d2i,  and  du  set  at  0.779,  0.949,  and 
0.880  A,  respectively,  a  12  was  varied  from  1.864  to  1.939 
A  in  steps  of  0.025  A  (see  runs  23—26  in  Table  II).  A 
parabolic  fit  to  the  r  factor  as  a  function  of  al2  gave  a 
minimum  for  a12  =  1.934  A.  The  next  series  of  calcula¬ 
tions  varied  dn,  d2},  and  di4  independently  in  small  steps 
near  the  above  rough  best  values  (0.779,  0.949,  and  0.880 
A,  respecively),  with  al2  fixed  at  1.934  A  and  a23  bulk¬ 
like.  Figure  4  shows  the  values  of  the  r  factor  for  varia¬ 
tions  of  d\2  and  d2}  set  to  0.880  A  (the  same  plot  could  be 
made  for  nearby  values  of  d 34,  but  the  r-factor  values  are 
slightly  higher).  Fitting  a  quadratic  function  of  dn  and 
d23  (elliptic  paraboloid)  to  the  r  factor  gives  dl2  =0.759 
A  and  d2i  =  1.014  A  as  optimum  values.  With  the  first 
three  interlayer  spacings  fixed  at  the  above  optimum 
values,  a  )2  was  varied  from  1.903  to  1.994  A  in  steps  of 
0.030  A  (runs  27—30  in  Table  II),  a  parabolic  fit  to  the  r 
factor  gave  ai2=  1.949  A.  Next,  with  the  above  best 
values  of  d\2,  d2i,  d^,  and  a12,  calculations  were  made 
for  a2j  changing  from  1.752  to  1.934  A  (runs  31—34  in 


TABLE  III.  Results  of  calculations  for  0 =  14*,  — 107.5’  of  various  structural  models:  d |2,  </2j, 

and  1/34  are  the  first  three  interlayer  spacings,  o)2  and  a2j  are  the  first  two  interlayer  registry  parame¬ 
ters  (see  text).  The  agreement  between  theory  and  experiment  is  given  by  the  reliability  factor  r.  All 
distances  are  in  A. 


Run 

dl2 

d2  j 

dM 

«if 

fl23 

r 

1 

..EOT 

0.869 

1.949 

EE 

2 

0.869 

1.949 

,.■81.. 

■Mia 

3 

0.758 

0.996 

0.869 

1.949 

1.853 

0.125 

4 

0.758 

1.012 

0.869 

2.002 

1.853 

0.137 

5 

0.758 

1.012 

0.869 

1.949 

l.&uO 

0.120 

6 

0.764 

1.018 

0.869 

1.936 

1.850 

0.120 

7 

0.770 

1.024 

0.863 

1.923 

1.846 

0.119 
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TABLE  IV.  Optimum  values  of  structural  parameters  and  of  the  real  part  of  the  inner  potential  K0. 
A£  equals  total  energy  range. 


Data 

dn 

da 

d* 

a  12 

023 

Vo 

A £  (eV) 

r 

0=0° 

0.759 

1.014 

0.880 

1.949 

1.841 

-10.5 

1704 

0.098 

<3> 

II  II 

i.5 

s 

LA 

« 

0.772 

1.026 

0.869 

1.917 

1.845 

-9.5 

1185 

0.119 

Mean 

0.764 

1.019 

0.874 

1.936 

1.843 

-10.1 

2889 

Table  II).  The  /--factor  minimum  occurs  for  a2j  =  1.841 
A.  The  best  value  of  the  real  part  of  the  inner  potential, 
Vo,  which  was  allowed  to  vary  independently  in  all  calcu¬ 
lations,  was  determined  to  be  — 10.5  eV. 

We  wish  to  emphasize  at  this  point  that  the  above 
method  for  finding  the  overall  r- factor  minimum  (and 
hence  the  solution  to  the  surface  structure  problem)  in  the 
six-dimensional  parameter  space,  with  r  =r<di2,d2), 
<^34>ai2>a23>*o)>  is  incomplete  owing  to  the  limited  num¬ 
ber  of  (lengthy)  calculations  which  were  feasible.  The 
main  problem  is  the  correlation  between  different  vari¬ 
ables,  i.e.,  it  is  not  possible  to  vary  one  parameter  at  a 
time  in  order  to  optimize  the  structure.  In  analyzing  the 
normal-incidence  data  for  Fe{310)  we  have  made  a  sub¬ 
stantial  number  of  calculations  allowing  independent  vari¬ 
ations  of  the  three  most  sensitive  parameters,  dn,  d2 3, 
and  c/34,  while  consideration  of  changes  in  a!2  and  a2 3 
were  more  limited.  In  order  to  take  into  account  more  ac¬ 
curately  the  correlation  between  registry  shifts  and  inter¬ 
layer  spacings,  a  different  strategy  was  used  for  analyzing 
the  data  at  0=14°,  #=  —  107.5°.  The  method  used  was 
the  following:  (i)  Use  values  of  the  structural  parameters 
set  near  the  optimum  values  found  from  analyzing  0=0° 
data  as  initial  values;  (ii)  vary  each  of  the  parameters  one 
at  a  time  by  a  small  amount;  (iii)  calculate  the  local  gra¬ 
dient  Vr  of  the  r  factor  as  a  function  of  the  structural 
variables;  and  (iv)  make  a  series  of  calculations  along  the 
direction  of  —  Vr  in  the  parameter  space  and  find  the  r- 
factor  minimum  along  that  direction.  If  the  surfaces  of 
equal  r  factor  are  spherical,  then  —  Vr  will  point  exactly 
toward  the  global  minimum.  Problems  arise  when  the 
equi- r-factor  surfaces  are  not  spherical  (and  they  usually 


TABLE  V.  r  factors  for  the  bulklike  and  fully  relaxed 
models  for  the  three  experimental  data  sets  and  mean  r  factors 
for  all  data. 


Data  sets 

Bulklike 

Relaxed 

Energy  range  (eV) 

b 

II 

0.292 

0.102 

1704 

(t-T 

#=-107.5° 

0.193 

0.126 

1205 

fl— 14° 

#=-107.5° 

0.187 

0.127 

1185 

Mean 

0.232 

0.116 

4094 

are  not)  and/or  when  the  r  factor  has  multiple  minima 
(this  case  is  also  common) — the  object  is  to  find  the 
deepest  minimum.  The  solution  to  the  former  case  is  ob¬ 
tained  by  recalculating  the  local  gradient  about  each  new 
minimum  until  the  point  is  found  where  the  local  gradient 
vanishes.  The  latter,  more  difficult,  case  involves  blanket¬ 
ing  the  parameter  space  with  enough  calculations  so  that 
the  region  with  lowest  absolute  values  of  the  r  factor  is 
found  and  then  local-gradient  method  can  be  applied.  It 
is  clear  that  when  there  are  many  variables,  compromises 
must  be  made  which  will  increase  the  uncertainties 
beyond  those  inherent  in  the  experimental  technique. 

Table  III  shows  the  calculations  used  to  determine  the 
local  gradient  of  the  r  factor  as  a  function  of  dl2  d2i,  a )2, 
and  a  23-  The  values  of  d^  and  V0  were  varied  indepen¬ 
dently  in  each  calculation  since  little  extra  computer  time 
was  required;  the  entries  in  the  table  are  for  the  optimum 
values  of  0.869  A  for  dH  and  —9.5  eV  for  V0.  Runs  1—5 
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FIG.  6.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  thr  fiilly 
relaxed  model. 
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energy  <«v> 


FIG.  7.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  the  fully 
relaxed  model. 


FIG.  9.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  the  fully 
relaxed  model. 
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ENERGY  («V) 

FIG.  8.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  the  fully 
relaxed  model. 


ENERGY  (»V) 

FIG.  10.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  the  fully 
relaxed  model. 
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ENERGY  («V) 

FIG.  11.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  the  fully 
relaxed  model. 


in  Table  III  were  used  to  evaluate  Vr  =  (di2,d23,fl|2,ay )> 
and  runs  6  and  7  were  made  along  the  direction  —  Vr. 
The  parameter  values  determined  by  a  parabolic  fit  to  the 
r  factors  of  runs  1,  6,  and  7  are  given  in  Table  IV  as  are 
the  values  from  the  0=0°  analysis  and  the  weighted  (by 
the  total  energy  range  of  the  data  sets)  mean  values.  The 
three  decimal  places  recorded  in  the  table  represent  the 
precision  of  the  computations,  not  the  accuracy  (see 
below).  Nevertheless,  the  agreement  between  the  two  data 
sets  is  very  good;  the  values  for  corresponding  interlayer 
spacings  are  within  0.01  A  of  each  other  and  the  registry 
shifts  are  consistent  within  0.03  A.  In  Table  V  we  list  the 
r-factor  values  for  the  data  at  three  angles  of  incidence 
for  the  bulklike  and  fully  relaxed  models.  The  spectra  at 
0=7”,  <f>—  — 107.5”  were  not  used  in  the  intensity  analysis 
owing  to  limitations  on  computer  time.  Figure  5  is  a 
schematic  representation  of  the  relaxed  surface  structure 
and  Figs.  6—12  show  the  experimental  spectra  together 
with  the  calculated  spectra  for  both  the  relaxed  and  bulk¬ 
like  structures-  Comparison  of  spectra  for  the  relaxed  and 
bulklike  models  with  each  other  and  with  experiment 
demonstrates  the  sensitivity  of  LEED;  the  overall  im¬ 
provement  in  agreement  between  theory  and  experiment 
obtained  for  the  relaxed  structure  is  quite  striking. 

In  summary,  the  observed  structure  of  Fe(  310}  is: 

(i)  bulk  interlayer  spacing  0.906  A, 

(ii)  bulk  layer-to-layer  registry  shift  1.813  A, 

(iii)  d|2=0. 76±0.3  A  (16. 1  ±3. 3  %  contraction). 


30  50  70  90  110  130  150  170  190 

ENERGY  (tv) 


FIG.  12.  Experimental  and  calculated  LEED  spectra. 
Theoretical  curves  are  shown  both  for  the  bulklike  and  t**e  fully 
relaxed  model. 


(iv)  d23  =  1.0210.03  A  (12.6±3.3  expansion), 

(v)  d34  =0.8710.04  AJ4.0±4.4%  contraction), 

(vi)  a,2  =  1.94+0.05  A  (7. 2 ±2. 8%  relaxation  toward 
more  symmetrical  registration  with  respect  to  the  second 
layer), 

(vii)  a23  =  1.84+0.05  A  (1.612.8 %  in  same  direction 
as  top-layer  registry  shift), 

(viii)  V0=  — 10.  lll.OeV, 

(ix)  rmjn  =0.116  (Zanazzi  and  Jona13  r  factor,  33  spec¬ 
tra  for  a  total  4094-eV  energy  range). 

The  percentage  changes  are  with  respect  to  the  bulk 
values  given  above.  The  associated  errors  are  estimates; 
for  a  discussion  of  errors  in  LEED,  see  Ref.  20. 

The  change  in  top-layer  #  registry  shift 
Aai2=a|2(relaxed)-a|2(bulk)=0.13  A,  is  to  be  compared 
to  the  change  of  0.23  A  which  would  place  each  first-layer 
atom  in  a  position  where  it  would  have  three  nearest 
neighbors  in  the  second  layer  as  opposed  to  two  for  the 
bulklike  registry.  Thus  the  parallel  shift  of  the  top  layer 
is  56%  of  the  shift  necessary  to  achieve  maximum  coordi¬ 
nation  of  first-layer  atoms  with  respect  to  the  second 
layer,  the  shift  being  limited  by  the  influence  of  deeper 
layers.  The  change  in  «econd  registry  shift  is  in  the  same 
direction  as  the  first,  in  contrast  to  the  case  of  Fe{211J 
(Ref.  10)  where  the  first  two  changes  in  registry  shift  were 
found  to  have  opposite  signs.  There  are  not  enough  cases 
known  at  present  of  parallel  relaxations  to  determine  the 
systematica. 
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